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Using the basis of Hermite-Fourier functions (i.e. the quantum oscillator eigenstates) and the 
Sturm theorem, we derive the practical constraints for a function and its Fourier transform to be 
both positive. We propose a constructive method based on the algebra of Hermite polynomials. 
Applications are extended to the 2-dimensional case (i.e. Fourier-Bessel transforms and the algebra 
of Laguerre polynomials) and to adding constraints on derivatives, such as monotonicity or convexity. 



I. INTRODUCTION 

Positivity conditions for a Fourier transform appear in various domains of Physics. Frequent questions are, 
• What are the constraints for a real function i/j(r) ensuring that its Fourier transform, 

& ^ (s) = ^T e is ^{r) , (1) 
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be real and positive? 

Conversely, what are the properties of ip if ip is positive? 



• Finally, what are the constraints on Fourier partners such that both ib and w be positive? 

^-H ' Physicists often need to work with concrete constructions. The practical construction of a basis of functions satisfying 
the abovementionned positivity properties remains, up to our knowledge, an open problem. Such questions are quite 
relevant in Physics. As practical examples, let us quote two well-known cases. A Fourier transform relates [1] two 
quantities, namely the cross section and the profile of a nucleus, which ought to be both positive. In Particle Physics, a 
2-d Fourier-Bessel transform relates the color dipole distribution in transverse position space (derived from Quantum 
Chromodynamics) and the transverse momentum distribution of gluons probed during a deep- inelastic collision [2]. 
In short, such questions occur also in probability calculus, for the relation between probability distributions and char- 
acteristic functions [3], in crystallography and more generally in condensed matter physics, e.g. for the interpretation 
' of patterns, etc.... 

The problem is simplified if related to another one, that concerning the functions which are invariant [4] up to a 
phase factor 1 by Fourier transforms. Indeed, the most familiar examples of positive self-dual functions or distributions, 
thus trivially verifying the double positivity condition, which are at most scaled under Fourier transformation (FT), 
are the Gaussian and the Dirac comb. Many special cases can be found, where positivity is conserved, such as, for 
instance, the continuous family of functions exp(— r v ), where < v < 2. Various sufficient conditions for positivity 
can be found in the literature, such as the convexity of ip [7] but, up to our knowledge, no general constructive method 
has been presented. 

The present note attempts to give general positivity criteria, in a constructive way, by taking advantage of a 
representation under which the FT is essentially "transparent" . Our method combines the advantages of self-duality 
properties with those allowed by an algebra of polynomials, where positivity means absence of real roots, hence 
explicit manipulations of polynomial coefficients. For this sake, in the 1-d case, we select a basis made of convenient 
eigenstates of the FT, the Hermite-Fourier functions, i.e. the harmonic oscillator eigenstates. The method extends to 
the 2-d case, or Fourier-Bessel transform, by replacing Hermite by Laguerre polynomials. 

There are general mathematical theorems about the characterization of Fourier transforms of positive functions [8] . 
Let us quote in the first place the Bochner theorem and its generalizations [9] which state that the Fourier transform 
of a positive function is positive- definite. But positive definiteness in the sense of such theorems does not imply plain 
positivity 2 . Hence our problem actually could be rephrased [10] as "build positive-definite functions that are positive". 
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1 They are called "self- Fourier" in [5] if the phase factor is 1 and "generalized self- Fourier" [6] (or "dual") for other phases. 
2 Positive definiteness means that for any real numbers xi...x k and complex numbers fi...£fc, one has . <p(x k —Xj) £j£ fe > 0. 
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Our formalism is the subject of Section II. Numerical, illustrative examples will be given in Section III. Then Section 
IV describes the extension of our algorithms to the 2-d problem. Hcrmite polynomials will be replaced by Laguerre 
ones, but the algebra remains essentially the same. A brief discussion, conclusion and outlook make Section V. 



II. BASIC FORMALISM 



Consider the harmonic oscillator Hamiltonian, \{p 2 + r 2 ), and its eigenwavefunctions, 



u n [r) = tt 4 e 2 



H„(r). 



(2) 



Here, we set H n to be a square normalized Hcrmite polynomial, with a positive coefficient for its highest power term. 
For the sake of clarity, we list the first polynomials as, H = 1, Hi = \/2 r, H 2 = (2r 2 - 1)/V2, H 3 = (2r 3 - 3r)/V3 
and their recursion relation 



(3) 



(4) 

and thus such states give generalized self-dual functions with phase i n . If one expands ip in the oscillator basis, 
VK 7 *) — J2n=o VVi u n(r), with a truncation at some degree N, then all odd components ip2 P +i must vanish if ip must 
be real, and the even rest splits, under FT, into an invariant part and a part with its sign reversed, namely 



a n +iH n+ i — 2ra n H n — 2na„_i_ff„_i , 
where a n — \Jl n n\ . It is known that the FT of such states brings only a phase, 

dr e tsr u n (r) = i n u n (s) . 



^( S )=7T-i e -5 s [P+(s) - P-(s)\, P+(s) 



[N/4] 

tp4 P H ip (s), P_ 

p=0 



[(JV-2)/4] 

( s ) = i>4p+2 H 4p+2 {s) , 

p=0 



(5) 



where the usual symbols [N/4] and [(N — 2)/4] mean, respectively, the entire parts of N/4 and (N — 2)/4. 

Notice that, when all components ipn vanish except ipo, then both ip and ip are positive, because Hq = 1. Hence, 
one may, starting from this special point in the functional space of functions, investigate those domains of parameters 
ip n where the polynomials V = P + + P- and Q = P + — P_ have no real root. Notice that only even powers of r and 
s are involved. It will therefore be convenient to use auxiliary variables such as p — r 2 and a — s 2 , and the domain 
of interest for the parameters ip n will correspond to the absence of real positive roots for both p and a. 

The second ingredient of our approach is the well-known Sturm theorem [11] which gives the efficient way to 
characterize and localize the real roots of any given polynomial. The Sturm criterion can be expressed as follows: 

" Given a polynomial V(x), its Sturm sequence S(x) = {<Si,<S2, ... S m , ...Sj<n} is the set of polynomials 



Si =V, S 2 = dV/dx, ... S m = -S m -2 + 



(6) 



where [-] designates the polynomial quotient 3 . To know the number of roots between x = a and x = b, count the 
number J\[(a) of sign changes in S(a) and, similarly, count N(b). Then the number of roots is \Af(b) — A/"(a)| . " 

The domain borders where the root number, \Af(+oo) — A/"(0)| , changes have to do with cancellations of the resultant 
1Z between V and dV/dx. The cancellation of 1Z corresponds to collisions between conjugate complex roots becoming 
real roots and conversely. Because of the demanded positivity of p and a, the borders have also to do with sign 
changes of V(0) and Q(0), meaning real roots p and a going through 0. All such technicalities are taken care of by 
the Sturm criterion, which furthermore allows the labeling of each domain by its precise number of real roots. 

This will be implemented here in an explicit way, analytically as much as possible, then numerically and graphically, 
for a few cases of a general illustrative value. For this, we will plot the shapes of domains labeled by the values of 
the Sturm criterion. Hence the combination of the self-dual properties of the quantum oscillator basis and of the 
Sturm theorem allows a constructive method for a systematic investigation of positivity conditions for a 1-d Fourier 
transform. This basis has the potential to represent any function in the Hilbcrt space C 2 , but the present study 
concerns a finite set of components. An extension to infinite series of components may deserve other tools. 



3 The Sturm sequence is thus made of polynomial remainders, with (— ) signs. It obviously stops at some j < N. 
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III. POSITIVITY DOMAINS 



In this Section we will apply our method to the case of a basis formed with 3 or 4 Hcrmite-Fourier functions. We 
consider: A the basis 0o, 04, ipg, then B the basis 0o, 04> 08j 0i2j which both lie in the subspacc with eigenvalue 1, and 
C the basis 0o,02, 04 where 02 is in the subspace with eigenvalue —1, furthermore, D the influence of an additional 
constraint motivated by Physics, that of monotony for ip, and finally E a comparison with the convexity constraint 
for 0. Many other illustrations are possible, but these, A-E, demonstrate the flexibility of our approach. 



A. Mixture of three polynomials in the subspace with eigenvalue 1 

Here we assume that has only components 0o, 04 and 0s, hence V reduces to P + in (5) and we can study 

v = ^ + 27E (4p2 " 12/9 + 3) + 24vW (16// " 224p3 + 840p2 " 840/3 + 105 } ' (7) 

where p = r 2 . One maintains 8 > and one reads -P+(0) = 0o + 304/(2 \/6) + 1050s/ (24 V70). The resultant 1Z 
between P and dP/dp is, 

K cx 2100 00 01 - 1050 V6 0^ - 240 V7O0^ 0^ 8 + 400 \/lO50o 0l 0s - 165 V7O0|0 8 + 
48O0^0| +4560 V6 0o 4 ~ 13320 O 04 0| - 16OO\/6 0f 0f - 792^70 00 0f + 
1728v / lO5 0o04 0f - 612v / 7O0l0| - 10080 o 0g - 252OV6 4 0g + 1260^01 . (8) 




FIG. 1. Mixture of Ho, Hi and H$. White, no real positive p root; black, 4 roots; darkening grays, increasing root number. 

It is convenient to take advantage of the free scaling of V and normalize the polynomial so that its coefficients lie 
on a half sphere of unit radius, 

0o = cos a cos (3, 04 = sin a cos (3, 08=sin/3, — n < a < ir, 0</3<7r/2. (9) 

We thus show in Figure 1 the domains where the number of roots increases from to 4. The no root domain is the 
white triangle above a = (3 = and slightly right from this point. In this domain, (= p) is self-Fourier and positive. 
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B. Mixture of four polynomials in the subspace with eigenvalue 1 



Now we add to ip a component ipxi > hence 

V = ib + (V - 12p + 3) + (16p 4 - 224p 3 + 840p 2 - 840p + 105 ) + 

2 V6 24 a/70 

' 12 : (G4p 6 - 2112p 5 + 23760p 4 - 110880p 3 + 207900p 2 - 124740p + 10395 ). (10) 



While borders corresponding to P+ (0) — obtain easily, the resultant 1Z to be considered for other borders is unwieldy 
and is skipped here. Taking advantage of scaling we set, 

ipQ = cos a cos (3 cos 7, ip4 = sin a cos/3 cos 7, Vs = sin /3 cos 7, V12 = sin 7, 

-|<a<|, -tt<(3<tt, 0< 7 <|. (11) 

This choice of spherical coordinates was designed to ensure the positivity of ipi2, obviously, but also a domi- 
nance of tjjQ near a = j3 = 7 = 0. The dominance is clearly useful for the positivity of V . Then this S3 sphere 
can be explored by various cuts according to fixed values of 7. The results are shown in Figures 2-5, with 
7 = 7r/10 6 , 7r/15, 27r/15, 7r/6, 7t/5, 77r/30, 47r/15, 37r/10, respectively. The color code for the number of roots is: 
root, red; 1, yellow; 2, yellowish green; 3, bluish green; 4, blue; 5, dark purple; 6, pink (in the uncolored edition, 
the no root domain, if it exists, is that dark, small or tiny triangle slightly right of the map center). The red domain 
shrinks at first very slowly when 7 increases, then faster when 7 ~ 7r/6. Beyond such an order of magnitude for 7, 
there is no red domain and the map becomes invaded by bigger and bigger pink patches, representing the dominance 
of the 6 positive, real roots of 



-6 

3.14159 10 




FIG. 2. Mixture of Ho, H4, H% and Hn- Root number maps. Color code: root, red; 1, yellow; 2, yellowish green; 3, 
bluish green; 4, blue; 5, dark purple; 6, pink. Left: cut of the parameter sphere S3 when 7 = tv/10 6 . Right: cut for 7 = 7r/15 
(uncolored edition: the root domain is the dark triangle-like domain near the center). 
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FIG. 3. Left: Same as Fig. 2, with 7 = 27r/15. Right: 7 = n/6. See how the no root triangle shrinks. 
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FIG. 4. Left: 7 = n/5; Right: 7 = 7n/30. Absence of no root domain. 
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FIG. 5. Left: 7 = 47r/15; Right: 7 = 37r/10. Progressive dominance of Hn with 6 p roots. 



C. Two polynomials from subspace "1" mixed with one polynomial from subspace "—1" 

If we consider a mixture of Hq, H2 and H±, the FT connects the two polynomials 

^, s , ,20-1 , V-12p + 3 x , , 2a -I , 4cr 2 -12cr + 3 

We study the positivity of each polynomial separately, then of both. Notice that the parametrization, 

-00 = cos a cos /3, tp2 — sin a cos /?, V4 = sm /3, — 7r < a < ir, 0</3<7r/2, (13) 

reverses only the sign of tp2 if ct becomes —a. This parity operation is seen in Figure 6, the white domains of which 
correspond to the positivity of V and Q, respectively. The domain of simultaneous positivity for both is the white 
intersection domain in the left part of Figure 7, with the expected symmetry. 

It is actually easy here to analyze analytically the resultants of interest for V and Q, 

ftoc Vo>2-4Vo^4 T 4V2> 2 ^4 + 2Vo>|, (14) 

together with signatures for the signs of roots, such as 'P(O), etc. This can be done also in the "spherical representa- 
tion". The white domains of Figs. 6, 7 are recovered. 



G 
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D. Positivity with monotony 



For some problems [2], it may be useful to request either ip and/or to be monotonous functions in an interval 
such as [0, oo [. We illustrate this in the case of an H , H 2 , H4 mixture, with the additional constraint, 

— oc — 12 cos a cos /3 + 6\/2cos/3sina(2cr-5)-sin/3V6(4cr 2 -28cr + 27) < . (15) 
da 

The result appears in the right part of Fig. 7. The white domain, corresponding to such three simultaneous conditions 
of positivity and monotonicity, is a severe restriction of the white domain seen in the left part of Fig. 7. 



E. Positivity from convexity 



A practical and sufficient, but not necessary condition for the positivity of ip is the convexity of ip [7]. Indeed, 
<p <~ Jq 00 dr [1 — cos(sr)]/s 2 d 2 ip/dr 2 > 0. We illustrate this convexity condition for a mixture of H , H4, H%. It is clear 
that the presence of exp(— r 2 /2) in front of a finite order polynomial, with the even parity of ip and its derivability, 
are contradictory with "convexity everywhere" ; a smooth, round maximum must occur at the origin. We thus study 
partial convexity conditions of the kind, "convexity between r c and +00". A reasonable choice for the order of 
magnitude of r c is the position r c = 1 of the inflexion point of exp(— r 2 /2) . The second derivative d 2 ip/dr 2 belongs to 
the same algebra. We adjusted its Sturm criterion to various values of r c . Small domains only are found that ensure 
zero roots for d 2 ip/dr 2 , because most mixtures of Ho, H4, H$ do oscillate. Figures 8 and 9 show, in white again, with 
the parametrization by Eqs. (9), the survivor domain obtained if r c = 1, right part Fig. 8, then r 2 = 2/3, left part 
Fig. 9, and r 2 = 2/5, right part Fig. 9, respectively. The domain shrinks in a smooth way when r 2 decreases from 1 
to 0.4 and disappears if r 2 <^ 0.4. It does not increase much when r c > 0(1). The left part of Fig. 8, a zoom of the 
left part of Fig. 6, shows that such partial convexity domains are already included in the positivity domain of ip. 




FIG. 8. Ho, H$, H& mixture. Left: zoom of the positivity domain. Right: convexity domain for r > r c = 1. 
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FIG. 9. Same mixture as Fig. 8. Convexity domains. Left: r\ — .67. Right: r\ = .4. 

IV. POSITIVITY FOR THE 2-DIMENSIONAL FOURIER TRANSFORM 

The Fourier-Bessel transform in which we are here interested reads, 

/>oo 

= / (r dr) J {sr) ip(r) . (16) 
Jo 

For the 2-d radial space, a complete basis of states results from substituting r 2 for r into Laguerre polynomials, 

roc 

/ (rdr) 2 e- r " L m (r 2 ) L n (r 2 ) = S mn . (17) 
Jo 

For the sake of clarity, we list here the first four such normalized, "2-d radial" states, 

_i r 2 f 2 r 4 -4r 2 + 2 r 6 - 9r 4 + 18r 2 - 6 1 
{v , v u v 2 , v 3 } = V2e * r il, r 2 -l, , |. (18) 

One can verify that the states v n make eigenstates of the Fourier-Bessel transform, 

/>oo 

/ (rdr)J (sr)v n (r) = (-l) n v n (s). (19) 
Jo 

Positivity conditions can again be implemented with the Sturm criterion. For instance a mixture of ^cb 1 ^, from that 
subspace with eigenvalue 1, and v\, from that with eigenvalue —1, defines the following two reciprocal partners, 

V{p)=1> Q + ih (P-1)+V>2 p2 ~t P + 2 , Q(<t)=Vo-^i {a-l) + ^~ A ° + 2 ■ (20) 




FIG. 10. Positivity domains for the "radial 2-d", parity mixed case. Left: ip. Right: both tjj and ip. 

The left part of Figure 10 shows the positivity domain for (p and its right part shows the joint positivity domain 
for ip and ip. Similarities beween polynomials involved in the present "mixed parity" case and those of case C of the 
previous Section create topological similarities with Fig. 6 (right) and Fig. 7 (left) but numerical details do differ. 




-3 -2 -1 Q 1 2 3 -0.1 -0.05 0.05 0.1 0.15 0.; 
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FIG. 11. Left: Positivity domain for the "radial 2-d", vq, V2, V4 mixture. Right: third derivative negativity domain if r\ = 1.4. 
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Interestingly enough, a sufficient condition [7] for the positivity of ip in this 2-d situation is that the third derivative, 
d 3 ip/dr 3 , be negative. But, as in the previous Section, case E, the presence of a truncated number of basis terms in 
our expansion may reduce this negativity condition for d 3 tp/dr 3 to a domain r > r t only. We show in Figure 11 the 
result for a mixture of vo, v<i and if r 2 = 1.4, a demanding situation resulting into the tiny domain in the right 
part of Fig. 11. The domain belongs to the positivity domain seen in the left part of the same Figure. 



V. DISCUSSION AND CONCLUSION 



To summarize our results, we have built a basis of functions verifying positivity together with their Fourier transform. 
The method is based upon algebras of Hermite polynomials (for 1-dimensional FT) or Laguerrc polynomials in the 
variable r 2 (for 2-dimensional radial FT). 

The Fourier transform has four eigenvalues l,i, — 1,— i, and thus four highly degenerate eigcnsubspaccs. Two of 
such subspaces are compatible with real functions remaining real. To span each subspace, we used a basis made of 
Hermite- Fourier (or "Laguerre- Fourier" ) states. The orders (in r) of the polynomials have to be multiples of 4 if the 
eigenvalue is 1 , and multiple of 4 plus 2 if the eigenvalue is — 1 . At the cost of a truncation of such bases to a maximum 
order N, the conditions for positivity, convexity, etc. thus reduce to simple manipulations of polynomial coefficients 
based on the Sturm theorem. In each truncation case, one can find suitable domains for the parameters which mix 
the various basis polynomials. Such domains have been illustrated by the figures shown by this paper. 




FIG. 12. Left: two examples of self Fourier states (ip = ip). Right: mixed parity case (ip / ip). 

Some qualitative considerations may be drawn about the solutions we found in various examples. In the left part 
of Figure 12, we display two typical solutions with self-Fourier properties. In the right part of Fig. 12 we display one 
solution connecting two distinct partners ip and <p. It is interesting to note that the shapes show distinctive features. 
One class of shapes, which are monotonic, seem to remain closer to the bare Gaussian. Indeed, the bare Gaussian 
is the building block of our method; it always belongs to the subsets we found. The other class, with oscillations, is 
different. One might be concerned with the possibility that our truncations to a maximum order N limit the flexibility 
of our method to the vicinity of the Gaussian. It seems not to be so, as shown for instance by the very oscillating 
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solution 4 in the left part of Fig. 12, which exhibits an approximate periodicity in an interval; it is reminiscent of the 
Dirac comb, which is of course outside the set of functions constructed with a finite number of polynomials. 

Our bases are flexible enough to reconstruct any function having positivity properties, but in some cases convergence 
might be slow. It is not excluded that other bases exist, that might be more convenient to speed up the convergence 
and make easier the search for positivity domains. Another open problem is that of positivity for periodic functions. 
Such questions are beyond the scope of the present paper. 
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